Dipole–dipole interactions enhance non-Markovianity and protect information against dissipation
Jan Munsif1, 2, Xu Xiao-Ye1, 2, Wang Qin-Qin1, 2, Chen Zhe1, 2, Han Yong-Jian1, 2, Li Chuan-Feng1, 2, †, Guo Guang-Can1, 2
Key Laboratory of Quantum Information of Chinese Academy of Sciences (CAS), University of Science and Technology of China, Hefei 230026, China
CAS Center for Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, China

 

† Corresponding author. E-mail: cfli@ustc.edu.cn

Project supported by the National Key Research and Development Program of China (Grant Nos. 2017YFA0304100 and 2016YFA0302700), the National Natural Science Foundation of China (Grant Nos. 61327901, 11474267, 11774335, and 61322506), the Key Research Program of Frontier Sciences, Chinese Academy of Sciences (Grant No. QYZDY-SSW-SLH003), the Fundamental Research Funds for the Central Universities, China (Grnat No. WK2470000026), the National Postdoctoral Program for Innovative Talents, China (Grant No. BX201600146), China Postdoctoral Science Foundation (Grant No. 2017M612073), and Anhui Initiative in Quantum Information Technologies, China (Grant No. AHY020100). The author Munsif Jan is thankful to the China Scholarship Council (CSC) for financial support (Grant No. 10358).

Abstract

Preserving non-Markovianity and quantum entanglement from decoherence effect is of theoretical and practical significance in the quantum information processing technologies. In this context, we study a system S that is initially correlated with an ancilla A, which interacts with the environment E via an amplitude damping channel. We also consider dipole-dipole interactions (DDIs) between the system and ancilla, which are responsible for strong correlations. We investigate the impact of DDIs and detuning on the non-Markovianity and information exchange in different environments. We show that DDIs are not only better than detuning at protecting the information (without destroying the memory effect) but also induce memory by causing a transition from Markovian to non-Markovian dynamics. In contrast, although detuning also protects the information, it causes a transition from non-Markovian to the Markovian dynamics. In addition, we demonstrate that the non-Markovianity grows with increasing DDI strength and diminishes with increasing detuning. We also show that the effects of negative detuning and DDIs can cancel out each other, causing a certain loss of coherence and information.

1. Introduction

Unwanted interaction between an open quantum system and its environment causes a detrimental effect called decoherence.[1,2] A variety of schemes have been proposed to overcome this problem, such as the quantum Zeno effect (QZE),[3] dynamical decoupling,[47] weak measurement and quantum measurement reversal.[8] However, previous theoretical studies have found that dipole–dipole interactions (DDIs) also help to reduce decoherence, in both two-atom[9] and three-atom[10] systems interacting with independent environments, protecting entanglement in both the Markovian and non-Markovian regimes. Non-Markovianity itself also helps to protect entanglement,[11] and similarly, various methods have been proposed for generating as well as preserving entanglement against decoherence.[1215] It has been found experimentally that, in the presence of DDIs, entanglement is sensitive to the number of trapped ions relative to the distance between them.[16] Zhang et al.[17] have reported that the DDIs show similar behavior to a positive detuning and exhibit the opposite role of negative detuning.

In cavity quantum electrodynamics,[18] the strong coupling regime (SCR), where the photons emitted by an atom can be reabsorbed and re-emitted, is a clear indication of information being exchanged between the system and environment,[19] implying that the process is non-Markovian (reversible process), and hence that the system can recover part of the lost coherence and information.[20,21] In recent decades, there has been increasing interest in the non-Markovian (reversible) nature of quantum processes,[2227] which has been found to be a vital resource in applications such as detecting the environment characteristic properties,[28] continuous-variable quantum key distribution,[2931] quantum metrology,[32] and steady-state entanglement.[33] Several studies have used the most common and significant measures of information flow to quantify the degree of non-Markovianity.[3439] To study the DDIs further experimentally, highly excited Rydberg atoms with large principal quantum numbers and strong DDIs are a good choice.[40] Currently, both Markovian and non-Markovian dynamics have been studied experimentally using ultra-cold Rydberg atoms,[41] and the transition from Markovian to non-Markovian dynamics has been realized experimentally in all-optical frameworks.[4244] In addition, Yuan et al.[45] have observed the transition from Markovian to non-Markovian dynamics caused by DDIs for a qubit coupled to a Bose–Einstein condensate (BEC) reservoir characterized by a pure dephasing spin-boson model.

In the recent schemes,[9,17] the authors limited their studies just to protect the entanglement via DDIs and detuning but in this study, we examine especially their impacts on the non-Markovianity (memory) and will see how they cause a transition from the Markovian to non-Markovian dynamics and vice versa. In detail, we consider a decoherence scheme,[46,47] where the system S is initially coupled to the ancilla A, interacting with the environment E via an amplitude-damping channel (ADC). We also go one step further and assume the realistic scenario where there are DDIs between the system S and ancilla A, while the system S and the environment E exchange information without direct interaction. Note that, during the processing of information, the survival of memory and entanglement plays an important role in the application of quantum technologies. In this regard, we have two parameters that can be used to adjust the system–environment coupling strength which controls the decoherence and dissipation. The first is the DDI strength, which increases the coupling between the system S and ancilla A, and the second is the detuning, which decreases the coupling between ancilla A and environment E.[48] Further, we analyze the role of DDIs and detuning, using the concepts of both accessible information [49] and entanglement of formation (EOF) to quantify the exchange of information between the system and environment. We show that entanglement sudden death (ESD) can be safely avoided and entanglement can be maintained for a long time without impacting the memory effect by properly tuning the DDIs in the non-Markovian regime. In our model, even weak DDIs ( ) are enough to cause a transition from Markovian to non-Markovian dynamics, and the degree of non-Markovianity (which we calculate numerically) is better than the case mentioned in Ref. [45]. In addition, comparing DDIs with detuning reveals that DDIs in our model provide better protection than detuning, and the proposals of any previous study.[9,10] We also analyze how negative detuning and DDI can cancel out the effects of each other. Ultimately, we hope that our model will be easier to realize in the laboratory than those of previous studies, due to its use of weaker DDIs.

The paper is organized as follows. In Section 2, we discuss the general behavior of Markovian and non-Markovian dynamics and also introduce the measures of non-Markovianity, such as accessible information , EOF . In Section 3, we present our example and its solutions in the presence of DDI. In Section 4, we present our results and discuss our findings. Finally, in Section 5 we present our conclusions.

2. Non-Markovianity measure via information flow

Generally, the quantum processʼs dynamical map is described by a time-local master equation[50] of the form

Here, the Lindblad superoperator is defined as
where H denotes the system Hamiltonian, ( ʼs), Fjʼs are the time-independent decay rates and Lindblad operators describing the effect of noise on the system, respectively. For positive ʼs, the dynamical map for Eq. (1) can be written as , satisfying the divisibility condition for
for all positive . This type of quantum dynamics is called as a conventional Markovian process. In time-dependent Markovian processes, the Hamiltonian H, noise operator Fj, and decay rates are explicitly time-dependent. If throughout the systemʼs evolution, then the dynamical map can be expressed in terms of a time-ordered exponential as , which transforms the state at time 0 to that at any other time t. Such dynamical maps have the fundamental property of satisfying divisibility condition. In such a scenario, the completely positive and trace-preserving (CPTP) map can be written as the composition of two other CPTP maps
where , for all . We also note that the value of may be temporarily negative during the systemʼs dynamics, in which the quantum dynamical map is indeterminate, which is no longer a CPTP map, and also violates the composition law for divisibility[51,52] given by Eq. (3).

2.1. Accessible information

Here, we discuss a well-known quantity that can extract the maximum amount of classical information about the system S by locally observing the environment E called accessible information,[49] is given by

where is interpreted as the von Neumann entropy of the reduced density operator of the system S and is a general projective measurement acting on the environmentʼs qubit E. The remaining state of the system S after measurement is given by , with probability in the subsystem E. Hence, any deviation of the accessible information from monotonically increasing behavior is a clear sign of non-Markovianity. Here, we demonstrate that, in the presence of DDIs, the dynamics remain non-Markovian for a long time and the system S exchanges some information with the environment E over time.

Non-Markovianity can be defined in terms of information flow via the accessible information. As we assume a zero-temperature reservoir, initially in a pure state, the composite tripartite state SAE always remain pure, and the Koashi–Winter relation[53] for the corresponding state can be characterized as

Here, represents the EOF shared by the system S and the ancilla A, which can be defined for the two-qubits case as
where
and
defines the concurrence,[54] in decreasing order of the eigenvalues , of the product matrix . In addition, , where is the complex conjugate of and is the Pauli spin operator. As we mentioned earlier, the system S and environment E do not interact directly, so the state is time-invariant throughout the dynamics after tracing over the subsystem A. Therefore, after taking the time derivative of Eq. (5), we obtain
This relation shows how information exchange occurs between the system S and environment E over time. It clearly indicates that any temporary increase in , during the open-system dynamics, leads to a temporary decrease in and vice versa. Thus, any temporary decrease in is a signal of divisibility violation and a direct consequence of non-Markovianity. However, for some non-divisible processes, may still decay monotonically. In that scenario, we can also use the entanglement-based measure of non-Markovianity criterion which can be expressed as
where the maximum is taken over all possible pure initial states of the bipartite system SA.

3. Example

Here, as shown in Fig. 1, we consider a system S that is initially entangled with an ancilla A, namely a two-level quantum system interacting with a zero temperature bosonic reservoir, represented by a collection of harmonic oscillators. We also consider the realistic scenario where there are DDIs between the system S and ancilla A, while the reservoir only influences the composite system SA via an ADC between A and E.

Fig. 1. A system S entangled with an ancilla A, and also including DDIs between them, while the ancilla A interacts with the environment E via an amplitude-damping channel (ADC).

The Hamiltonian that describes the entire composite system SAE is given by

where and denote the qubits raising and lowering operators, respectively, with transition frequency of . The environment annihilation (creation) operators are represented by bk , respectively, with mode frequencies . For continuous distribution of reservoir modes, the discrete coupling constant must be , where is the Lorentzian spectral density. Here, is the detuning of and the cavity central frequency, while λ defines the spectral width of the coupling; Both λ and are connected to the reservoirʼs correlation time and the systemʼs relaxation time via and , respectively. The static part of dipole–dipole interactions between system S and ancilla A is proportional to , where defines the relative position and with are the electric dipole moment of their respective dipoles. The weak-coupling regime, where the dynamics could be Markovian, corresponds to . In our scheme, the bipartite state SA is initially in a pure state, and the environment only interacts with the ancilla A, to gain information about the system S. Furthermore, we assume that the system S and ancilla A are maximally-entangled, interacting with a vacuum zero temperature bosonic reservoir which can be expressed as
The interaction between A and E induce ADC, that can be described by a unitary operator U as, and . In the operator sum representation the state dynamics can be formulated as, , where represent the corresponding time-dependent Kraus operator satisfying the relation for all time t,
Here, p(t) is the damping parameter which satisfies the relation
With the help of reservoir correlation function
and pseudomode approach,[9] we can easily calculate the damping parameter p(t) as
where is the inverse Laplace transform and η are the strength of DDI (detuning), respectively. In our case the time-dependent decay rate can be calculated as . The whole density state can be written as
Here, H.c represents the Hermitian conjugate of the last three terms. Now we trace out the environment E and the ancilla A, one by one, respectively, their corresponding reduced density states take the form
where , and , respectively.

4. Results and discussion

In this section, we study the effect of DDIs on the correlations dynamics in terms of EOF and accessible information as quantifiers. In the absence of DDIs, Figure 2(a) demonstrates the dynamics in the non-Markovian regime,[38] where we can observe that whenever increases temporarily, a decrease in leads to a negative decay rate , as shown in the inset of Fig. 2(a). In this case, the dynamical map is non-divisible, indicating a non-Markovian (reversible) quantum process. Next, we show that introducing DDIs between the system S and ancilla A, not only avoid ESD in the non-Markovian regime but also protect the information for a very long time. In Figs. 2(b) and 2(c) the DDI strengths are η=0.5, 1, respectively, which is still very weak relative to the strength, as reported in Refs. [9] and [10]. In our case, these plots show that information is exchanged between the system S and environment E much more rapidly. We can also observe (insets of Figs. 2(b) and 2(c)) that oscillates in the same way throughout the whole dynamics and is periodically negative, implying the quantum processes are non-divisible.

Fig. 2. Accessible information (blue line), the entanglement of formation (red line), for (a) η=0, (b) η=0.5, (c) η=1, respectively, with Δ=0, and in the insets, the decay rate (orange line) as a function of the dimensionless quantity , in the non-Markovian regime where . These plots demonstrate how information can be protected for a very long time without affecting non-Markovianity.

Similarly, Figure 3(a) represents the dynamics for zero DDIs in Markovian regime,[38] where we can notice that decays monotonically in this case, increasing and leading to positive throughout the evolution, as shown in the inset of Fig. 3(a). As long as remains positive, the dynamical map is divisible, indicating a Markovian quantum process. If we introduce DDIs in this Markovian regime, we observed that they disturb this monotonic behavior, depending on the DDI strength: Figure 3(b) shows that the behavior of and is no longer monotonic. This change, and the oscillating behavior of which is periodically negative throughout the dynamics, as shown in the inset of Fig. 3(b) is a clear indication of memory. It is important to note that, initially Markovian environment shows non-Markovian dynamics with the introduction of weak DDIs. We also demonstrate that the degree of non-Markovianity directly depends on the DDI strength. Increasing this causes information to be exchanged between the system S and environment E much more rapidly, protecting the information for a little longer, as we can see in Fig. 3(c). It is thus clear that, increasing the DDI strength in the strong-coupling regime causes the system and environment to exchange information faster without affecting its memory, whereas introducing weak DDIs causes a transition from Markovian (memory-less) to non-Markovian (with memory) dynamics. This is one of the main consequences of DDIs in the Markovian regime.

Fig. 3. Accessible information (blue line), the entanglement of formation (red line), for (a) η = 0, (b) η = 1, (c) η = 3, respectively, with Δ = 0, and in the insets, the decay rate (orange line) as a function of the dimensionless quantity , in the Markovian regime where . Here, panel (a) shows the monotonicity of information flow in the absence of DDIs, while panels (b) and (c) show non-monotonic behavior, revealing how the dynamics change when DDIs are introduced in the Markovian regime.

We also investigate the DDIs ability to protect entanglement, measured in terms of the concurrence. Recently, the authors in Refs. [9] and [10] have investigated using strong DDIs to prevent decoherence, whereas we have found that weak DDIs are sufficient to reduce entanglement degradation in our scheme. Figure 4 compares the behaviors of both DDIs and detuning in the non-Markovian regime, for different values of η (Fig. 4(a)) and Δ(Fig. 4(b)). From this we can conclude that DDIs are better than detuning at protecting entanglement in both regimes, but stronger DDIs are needed for protection in the Markovian regime than in the non-Markovian one. In addition, Figure 4(a) shows that DDIs can protect entanglement without destroying the memory effect. We can observe oscillations here for all DDI strengths, indicating that information is being exchanged between the system and environment, as also shown in the inset of Figs. 2(b) and 2(c). In contrast, protection through detuning leads to a loss of coherence, which in turn causes a transition from non-Markovian to Markovian dynamics. Although both DDIs and detuning provide less protection in the Markovian regime compared with the non-Markovian one, the important point here is the oscillation present in the DDI case, which represents a transition from Markovian to non-Markovian dynamics. Based on our present model, we prefer DDIs over detuning for two reasons: they provide better protection, preserving non-Markovianity in the strong-coupling (non-Markovian) regime, and induces memory in the weak-coupling (Markovian) case.

Fig. 4. Concurrence as a function of dimensionless quantity for (a) DDI values of η=0, 0.3, 0.6, 1, with Δ=0 and (b) detuning Δ=0, 0.3, 0.6, 1, with η=0, corresponding to solid red, black, blue, cyan lines, respectively, in the non-Markovian regime where . This comparison demonstrates that DDIs are better than detuning at protecting entanglement against dissipation.

Non-Markovianity enables information to return from the environment to the system. Figure 5 compares its behavior in the presence of DDIs and detuning. DDIs increase the coupling between the system and ancilla, whereas detuning reduces the coupling between the ancilla and environment, so although they both have the ability to protect information against dissipation but in case of non-Markovianity both have opposite behavior. Here, we use Eq. (9) to calculate numerically the degree of non-Markovianity for both DDIs and detuning over the whole time evolution. This shows that the non-Markovianity in the non-Markovian regime increases with DDI strength and decreases with detuning, corresponding to the red and blue dotted curves, respectively, in Fig. 5(a). Initially, the non-Markovianity is zero in the Markovian regime (Fig. 5(b)) but it arises with increasing the strength of DDI, representing a transition from a divisible to a non-divisible quantum map. We are able to achieve a much higher degree of non-Markovianity during this transition than as showed in Ref. [45]. From this, it is clear that increasing the coupling between the system and ancilla can protect the information and increase non-Markovianity while decreasing the coupling between the ancilla and environment can also protect the information but degrade the non-Markovianity.

Fig. 5. Non-Markovianity as a function of DDI (η) (red curve) and detuning ( ) (blue curve) throughout the whole time evolution for (a) non-Markovian regime with and (b) Markovian regime with . These demonstrate that non-Markovianity increases with DDI strength in both regimes, and decreases with increasing detuning in the non-Markovian regime.

As we have discussed earlier, DDIs show behavior opposite to that of negative detuning, although both are internal effects and can be used to prevent decoherence. Next, we investigate how negative detuning can be used to cancel out the effect of DDIs in the non-Markovian regime. Figure 6(a) shows that adjusting the strengths of both DDI and negative detuning to be equal but opposite can cause their effects to cancel out. We can also observe some non-periodicity for both weak DDIs and negative detuning, which can be eliminated by increasing their strengths in the same way. If we compare Fig. 6(a) ( ) with Fig. 2(a) ( ), we can see that although the figures are not identical, they do show the same behavior. Detuning causes a loss of non-Markovianity while protecting the information, whereas DDIs maintain it. Similarly, when DDIs and negative detuning are combined, they can cancel out each other, some coherence and information are lost, as can be seen from the amplitudes and numbers of peaks shown in Fig. 6(a). Further, intersection points of curves along time-axis clarify that these effects cannot cancel each other completely. In the non-Markovian dynamics, which possesses memory mean, future states depend on past states, introducing a weak detuning can eliminate this behavior. In Fig. 6(b) (η=0, Δ=0.5), there is an initial oscillatory phase, representing non-Markovian behavior, which later on reduces to monotonic behavior, indicating Markovian dynamics. The inset of Fig. 6(b) also reveals that there is a transition at , from a non-divisible to a divisible quantum process. The transition time also depends on the detuning strength. In summary, we conclude that DDIs cause a transition from Markovian to non-Markovian dynamics, while detuning can cause a slow and gradual shift from non-Markovian to Markovian dynamics, depending on the strength. We hope that this DDI trait will play a role in the quantum electrodynamics of atom–cavity systems against dissipation.

Fig. 6. Accessible information (blue line), the entanglement of formation (red line) for (a) , (b) η=0, Δ=0.5, respectively, and in the insets, the decay rate (orange line) as a function of dimensionless quantity in the non-Markovian regime where . Here, panel (a) demonstrates that the effects of DDIs and negative detuning can cancel each other out, while panel (b) illustrates the transition from non-Markovian to Markovian dynamics caused by detuning while it protects the information.
5. Conclusion

In this article, we have analyzed the decoherence scheme in the presence of DDI and detuning, in both Markovian and non-Markovian environments. Especially, we have considered a system S that is initially coupled to an ancilla A, where the environment E can gain information about the system S through ancilla A, without disturbing the system S directly. We have shown that, in the non-Markovian regime, even weak DDIs are sufficient to avoid ESD and enable entanglement to be maintained for a very long time. In order to prevent entanglement degradation, we need strong DDIs in the Markovian regime, but weak DDIs are sufficient to cause a transition from Markovian to non-Markovian dynamics. We have also found that the degree of non-Markovianity increases with DDI strength, and decreases with increasing detuning. It is also important to note that DDIs are more effective than detuning, both in protecting against decoherence and producing memory in the Markovian regime, whereas detuning causes a loss of memory (while still protecting the information) in the non-Markovian regime. In addition, we have found that the effects of DDIs and negative detuning can cancel out each other, with a certain loss of coherence and information. We believe that these results will be useful in understanding the role of DDIs and detuning to prevent decoherence, and prove to be helpful when realizing many future quantum protocols in practice.

Reference
[1] Lidar D A Chuang I L Whaley K B 1998 Phys. Rev. Lett. 81 2594
[2] Shor P W 1995 Phys. Rev. A 52 R2493
[3] Maniscalco S Francica F Zaffino R L Lo Gullo N Plastina F 2008 Phys. Rev. Lett. 100 090503
[4] Du J Rong X Zhao N Wang Y Yang J Liu R B 2009 Nature 461 1265
[5] Singh H Arvind Dorai K 2018 Phys. Rev. A 97 022302
[6] Peng S J Xu X K Xu K B Huang P Wang P F Kong X Rong X Shi F Z Duan C K Du J F 2018 Sci. Bull. 63 336
[7] Gong B Tu T Zhou Z Q Zhu X Y Li C F Guo G C 2017 Sci. Rep. 7 18030
[8] Wang S C Yu Z W Zou W J Wang X B 2014 Phys. Rev. A 89 022318
[9] Li Y Zhou J Guo H 2009 Phys. Rev. A 79 012309
[10] Khan S Jan M 2016 Int. J. Theor. Phys. 55 1515
[11] Rosario A Massoni E De Zela F 2012 J. Phys. B: At. Mol. Opt. Phys. 45 095501
[12] Huang L Y Fang M F 2010 Chin. Phys. B 19 090318
[13] Liao X P Fang M F Fang J S Zhu Q Q 2014 Chin. Phys. B 23 020304
[14] Yang G Lian B W Nie M 2016 Chin. Phys. B 25 080310
[15] Golkar S Tavassoly M K 2018 Chin. Phys. B 27 040303
[16] Yazdanpanah N Tavassoly M K 2016 J. Opt. Soc. Am. B 33 382
[17] Zhang Y Q Tan L Barker P 2014 Phys. Rev. A 89 043838
[18] Miller R Northup T E Birnbaum K M Boca A Boozer A D Kimble H J 2005 J. Phys. B: At. Mol. Opt. Phys. 38 S551
[19] Thompson R J Rempe G Kimble H J 1992 Phys. Rev. Lett. 68 1132
[20] Bellomo B Lo Franco R Compagno G 2007 Phys. Rev. Lett. 99 160502
[21] Franco R L Bellomo B Maniscalco S Compagno G 2013 Int. J. Mod. Phys. B 27 1345053
[22] Piilo J Maniscalco S Härkönen K Suominen K A 2008 Phys. Rev. Lett. 100 180402
[23] Tang J S Li C F Li Y L Zou X B Guo G C Breuer H P Laine E M Piilo J 2012 Europhys. Lett. 97 10002
[24] Xu J S Sun K Li C F Xu X Y Guo G C Andersson E Franco R L Compagno G 2013 Nat. Commun. 4 2851
[25] Fanchini F F Karpat G Castelano L K Rossatto D Z 2013 Phys. Rev. A 88 012105
[26] Chruściński D Maniscalco S 2014 Phys. Rev. Lett. 112 120404
[27] Addis C Bylicka B Chruściński D Maniscalco S 2014 Phys. Rev. A 90 052103
[28] Smirne A Cialdi S Anelli G Paris M G A Vacchini B 2013 Phys. Rev. A 88 012108
[29] Vasile R Olivares S Paris M A Maniscalco S 2011 Phys. Rev. A 83 042321
[30] Tan Y G Liu Q 2016 Chin. Phys. Lett. 33 090303
[31] Zhang S L Jin C H Guo J S Shi J H Zou X B Guo G C 2016 Chin. Phys. Lett. 33 120302
[32] Chin A W Huelga S F Plenio M B 2012 Phys. Rev. Lett. 109 233601
[33] Huelga S F Rivas A Plenio M B 2012 Phys. Rev. Lett. 108 160402
[34] Wolf M M Eisert J Cubitt T S Cirac J I 2008 Phys. Rev. Lett. 101 150402
[35] Breuer H P Laine E M Piilo J 2009 Phys. Rev. Lett. 103 210401
[36] Rivas A Huelga S F Plenio M B 2010 Phys. Rev. Lett. 105 050403
[37] Luo S Fu S Song H 2012 Phys. Rev. A 86 044101
[38] Fanchini F F Karpat G Çakmak B Castelano L K Aguilar G H Farías O J Walborn S P Ribeiro P H S de Oliveira M C 2014 Phys. Rev. Lett. 112 210402
[39] Haseli S Karpat G Salimi S Khorashad A S Fanchini F F Çakmak B Aguilar G H Walborn S P Ribeiro P H S 2014 Phys. Rev. A 90 052118
[40] Ravets S Labuhn H Barredo D Béguin L Lahaye T Browaeys A 2014 Nat. Phys. 10 914
[41] Genkin M Schönleber D Wüster S Eisfeld A 2016 J. Phys. B: At. Mol. Opt. Phys. 49 134001
[42] Liu B H Li L Huang Y F Li C F Guo G C Laine E M Breuer H P Piilo J 2011 Nat. Phys. 7 931
[43] Smirne A Brivio D Cialdi S Vacchini B Paris M G A 2011 Phys. Rev. A 84 032112
[44] Bernardes N K Cuevas A Orieux A Monken C Mataloni P Sciarrino F Santos M F 2015 Sci. Rep. 5 17520
[45] Yuan J B Xing H J Kuang L M Yi S 2017 Phys. Rev. A 95 033610
[46] Zurek W H 1981 Phys. Rev. D 24 1516
[47] Zurek W H 2003 Rev. Mod. Phys. 75 715
[48] Li J G Zou J Shao B 2010 Phys. Rev. A 81 062124
[49] Henderson L Vedral V 2001 J. Phys. A: Math. Gen. 34 6899
[50] Gorini V Kossakowski A Sudarshan E C G 1976 J. Math. Phys. 17 821
[51] Laine E M Piilo J Breuer H P 2010 Phys. Rev. A 81 062115
[52] Breuer H P 2012 J. Phys. B: At. Mol. Opt. Phys. 45 154001
[53] Koashi M Winter A 2004 Phys. Rev. A 69 022309
[54] Wootters W K 1998 Phys. Rev. Lett. 80 2245